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Abstract. In these proceedings we present a mini-review on the topic of the Dyson- 
Schwinger/Bethe-Salpeter approach to the study of relativistic bound-states in physics. In 
particular, we present a self-contained discussion of their derivation, as well as their truncation 
such that important symmetries are maintained. 


1. Introduction 

In Quantum Chromodynamics (QCD) the only observable objects are hadrons, which appear 
as bound-states of the elementary (but not observed) quark and gluon degrees of freedom. 
Consequently, most phenomenological aspects of QCD are essentially non-perturbative problems 
and require an appropriate framework for their study. One such approach is the combination 
of Dyson-Schwinger (DSE) and Bethe-Salpeter (BSE) equations that provide a means to study 
the non-perturbative properties of hadrons - at the microscopic level - without abandoning 
a priori the principles of QCD as a scale dependent continuum quantum field theory. We 
devote these proceedings to a concise exposition of known results about how this framework can 
be systematically constructed as well as a description of the main technical issues faced upon 
solving the DSEs and BSEs in combination. 

We refrain, for reasons of space, from including here results for observable data that can 
already be found in the literature [1-3]. Meson spectra have been thoroughly studied in [4-28] 
and references therein. Baryons have been studied in the quark-diquark approximation [29-37], 
and also as three-body objects [38-43]. Other observables of enormous experimental interest in 
deep inelastic scattering such as PDAs, PDFs, GPDs are also starting to be investigated [44-48]. 
Most of these calculations were performed upon truncating the (anti)quark-quark interaction 
to a single gluon exchange, the so-called Rainbow-Ladder truncation. This can be seen as the 
leading order interaction mechanism in the systematic expansion to be described below. While it 
provides accurate results for ground-state pseudoscalar and vector mesons, and also for ground 
state baryons, it shows clear deficiencies in all other channels. For this reason, current interest 
lies is the inclusion of interaction mechanisms beyond the leading one [17,49-63], and in the 
extension to glueball [64-66] and tetraquark bound-states [67] . 

The main purpose, therefore, of this contribution is to make it apparent that one is not 
pursuing a blind hunt for missing interaction mechanisms, but that there is a step-by-step 
programme for their inclusion in a systematic manner. 



2. Dyson—Schwinger and Bethe—Salpeter equations 

A complete description of a continuum quantum field theory, and in particular of QCD, is 
given when all the (infinitely many) Green’s functions of the theory are known. Functional 
methods deal with these without abandoning the realm of continuum physics, thus providing 
an approach complementary to that of lattice calculations. In particular, Dyson-Schwinger 
equations constitute an infinite set of coupled, non-linear integral equations for the full Green’s 
functions of the theory. We define in this section the basic concepts to be used later on. 

2.1. Functional Methods 

The generating functional for full Green’s functions corresponding to a classical Euclidean action 
S[(f>] is given by 

Z[J\ = f = (exp (—iJi4>i)) , (1) 

where denotes f d 4 xJ a (x)(j) a (x ) with i a superindex that subsumes possible discrete 

(collectively denoted a) and continuous indices x. This functional is normalised such that 
Z[ 0] = 1. The (•) indicate the weighted functional average with sources subsequently set to 
zero. In the presence of external sources one writes (-)j. 

The generating functional for connected Green’s functions is given by 

W[J] = —i\nZ[J] , (2) 

which can be shown by a Taylor expansion. Furthermore, the generating functional for one- 
particle irreducible vertex functions (the effective action) is obtained by a Legendre transform 

r[0 c ] = w[j\ - tfJi , with & = <&> j = 5 FF1 . (3) 

2.2. Dyson-Schwinger Equation 

The observation that the integral of a total derivative vanishes, true so long as the functional 
measure is invariant under spacetime translations of the field variables 

(4) 

can be used to obtain the Dyson-Schwinger equations (DSEs). Equivalently, one can view the 
DSEs as a consequence of the Ward-Identity associated with translational invariance. These 
DSEs are the quantum field theoretic equivalent of the classical equations of motion. Upon a 
vertex expansion we generate the infinite tower of non-linear integral equations that relate the 
fundamental Green’s functions of the quantum field theory to one-another. In a very dense 
notation [68] the DSEs for proper n-point Green’s functions may be obtained from 

sr[<i>] ss r 5 2 w 5 

~ Wi i + SJ5J54) 
by taking vacuum expectation values (•) of the nth derivative of the functional. 

2.3. Bound state equations 

Consider an n-particle Green’s function G^ n \pi ,..., p n ) describing the evolution of an n-particle 
system (each carrying momentum pi), and its amputated counterpart the scattering scattering 
matrix T^(pi,... ,p n ) 



G (n) = G (n) + G (n) r ( n ) G (n) , 


(6) 







with the disconnected product of n-propagators. Then we may obtain T from the Dyson 
equation 


T ( n ) =K + KG^T^ , 


( 7 ) 


with K the 2, 3,-, ..., n-particle irreducible interaction kernel. When the system forms a bound 
state, the momentum-space function T^ develops a pole at P 2 = —M 2 (in Euclidean spacetime) 
with P p = Pi the total momentum. At the pole we define 


rjn(n) 






P 2 + M 2 ’ 


( 8 ) 


where M is a state-dependent normalisation factor, is the bound-state Bethe-Salpeter 
amplitude and # its charge conjugate. Inserting this ansatz into the Dyson equation for the 
T-matrix and equating residues yields the homogeneous Bethe-Salpeter equation 


V = KG^V, (9) 

or for its conjugate 

^R- 1 = . (10) 


It can also be shown that the correct normalisation condition for Bethe-Salpeter amplitudes 
is [69,70] 




dG^ ] 

dP 2 



= 1 . 


( 11 ) 


Finally, Bethe-Salpeter amplitudes are the amputated versions of the Bethe-Salpeter wave 
functions <p = Gq^T. Using this notation, we can rewrite (9) as 


V = K<p, (12) 

The interaction kernels K contain a sum of infinitely many terms. If they are known, one can 
then solve the BSEs and obtain all the information about the bound states of the theory. 

3. Deriving the Bethe-Salpeter kernel 

The effective action provides a systematic means to derive (generalised) Bethe-Salpeter kernels. 
This not only enables new interaction mechanisms to be included in a controlled way, but ensures 
that relevant symmetries are maintained even when approximations are made. Here, we outline 
the main ideas behind this procedure for the case of two- and three-body bound-state equations 
following Refs. [71-73]. 

3.1. Local case 

We start with the case of local fields, introducing the main ideas applicable to the bi- and tri¬ 
local fields necessary for the study of bound-state equations. First, let us consider the generating 
functional for a purely fermionic action and introduce source terms Jk(x) for a local operator 
Ok(x) — {^(x), ^(x) 2 ,...}, with k an index describing all discrete indices of the operator 

and, also, distinguishing the different possible local operators 

Z[J] =M j Jk{ - x) ° k{x) ) , 


(13) 





where M = l/Z[J]\j_ 0 . Defining, as above, the generating functional for connected Green’s 
functions as W[J] = —ilog Z[J], the one-particle irreducible (1PI) effective action is given by 
the Legendre transformation 


r[*] = w[j\ - J d 4 xj k (x)* k (x ), 


(14) 


with the classical (super)field defined as 'I'j.(.x) = 5W[J]/SJ k (x). It can be determined in the 
absence of sources through satisfaction of the stationary condition 


<5r['h] 


= o. 


Consider a solution of this stationary condition, \l/£ 0 ^ (m). We call it stable if there exists solution 


(15) 


< 1} (*) = + A 'P k (x ) with A'lb (x) suitably infinitesimal. Expanding around (x 

yields 

6T[V] 


,(o) 


(o). 


Mb 


sr[v] 






( 0 ) 


— / d 4 x' 






6* k (x)6* e (x') 


*k=*l 0) 


A^(x') + • • • = 0 , (16) 


with a minus sign from exchanging the order of derivatives. To first order this equation implies 


J d 4 x' 


6* k (x)6*e (x f ) 


*k=*l 0) 


Affix') = 0 . 


(17) 


That is, the solution \H^(x) is stable if A^ k (x) is an eigenvector of r) 2 T/fhM'I' with zero 
eigenvalue. To interpret this, let us recall the following identity 


/ 


dfix" 


<5 2 r[^] 


8 2 W[J] 


■77 = -<W<5 (4) (x - x') 


8V e (x)8V k (x") 6J k (x")8J m (x') 

Now, consider an eigenvector £ k (x) of 8 2 W/8J8J with eigenvalue 1/A 

fi‘ X ' S2WlJl 

J x 6J k (*)6Jt{ X T X ’ 

Using the identity (18) we arrive at 


(18) 


(19) 


( 20 ) 


Therefore, £k(%) is also an eigenvector of 5 2 T/S^S^ with eigenvalue A. In particular, the 
perturbation AT^(x) corresponds to an eigenvector with A = 0, wich in turn is related to 
a pole in 5 2 W/5JSJ. In summary, stable solutions for are associated to poles of the 

connected Green’s functions given by 5 2 W/SJSJ (e.g. propagators if 0{x) = x )). 


3.2. Bilocal case: Two-body bound-states 

For the study of two body problems we add a bilocal term K^ s (x,y)0^ s (x,y) to the partition 
function (13), with k again distinguishing the different bilocal operators and the other discrete 
indices described by r, s. Specialising to the case of fermion-antifermion bound states (mesons) 
we limit the discussion to O rs (x,y ) = r ip r (x) r ip s (y) and antisymmetric sources K rs (x,y) = 














—K sr (y,x). Then, in addition to the classical field ^k( x ) we have the connected two-point 
Green’s function 

IggL = (21) 

The 2PI effective action is then defined as 

G] = W[j, K}- J d 4 xJ k {x)^ k {x) - t J d 4 xd 4 yK ab (x , y)tf 0 (®)tf 6 (y) (22) 

~\J d 4 xd 4 yK ab (x, y)G ab {x, y) , (23) 

which admits, in addition to (15) and in the absence of sources, the stationary condition 

■fMUo. (24) 

8G ab (x,y ) 


This is equivalent to the Dyson-Schwinger equation for the fermion propagator, as we show 
below. 

Following the ideas of previous section, a solution G^(x,y) of (24) is stable if there exists a 
perturbed solution G^(x,y) = G^}(x,y) + A G ab (x,y) determined by a non-trivial solution of 


d A x r d A y r 


5 2 T[T,G] 


8G ab (x,y)5G a 'v(x',y') 


G=G(°) 


AG a r b r(x' ,y') = 0 


(25) 


We show below that this is indeed equivalent to the usual Bethe-Salpeter equation for a fermion- 
antifermion bound state. Nevertheless this can be seen here, analogous to the local case, through 
the fact that a solution of (25) is related to a pole in S 2 W/5K5K , viz. to a pole in a four-point 
Green’s function. 


3.2.1. Relation to the Bethe-Salpeter equation It is known [74] that the 2PI effective action can 
be written as 


r[$,G] = S[tf] +iTr log G-iTrGo 1 G + r 2 [*,G\ , ( 26 ) 


where T 2 [T, G] contains two-particle irreducible diagrams only and Go is the classical propagator. 
The stationary condition for G gives, as promised above, the Dyson-Schwinger equation for the 
propagator. Indeed, taking a functional derivative with respect to the propagator G 


<sr[^,G] 

SG ab (x,y) 


iG ab(. X ^y)- iG 0^ab(. X ^y) + 


<5r 2 [fr,G] 

SG ab (x,y) 


= 0, 


(27) 


where we used 8Tr\ogG = G l 8G. Defining the self-energy as E = —idTo/dG we can rewrite 
the stationary condition for G as 


G db( x >y) = G oib( x >y) ~ ^ab(x,y) , 


(28) 


which is the gap equation for the fermion propagator. If G^ 0 ) is one of the solutions of the gap 
equation, we can take one further functional derivative with respect to the propagator G and 
rewrite (25) as 


= t d 4 x'd 4 y' — - - A G a ' b '{x',y') 

J 8G ab (x,y)8G a i b i{x',y') G=G0 

= J d 4 x'd 4 y' (-G^ (x , x')G$ (y, y') + K ab . cd (x, y; x', y')) A G cd {x', y'), 


(29) 










where we used SM^ 1 /SM^i = This is precisely the Bethe-Salpeter equation (12) 

for a Bethe-Salpeter wave function AG, with the quark-antiquark interaction kernel given by 


K a b-cd(xiy,x',y') = 


SG cd (x',y') 


G=G(°) 


(30) 


One then sees immediately that the interaction kernel is obtained by functionally cutting 
propagator lines from the self-energy. Note that the solution G = G^ is inserted into the 
self-energy only after the cutting has been performed. 


3.2.2. Goldstone bosons. Chiral symmetry. We show here that when the meson BSE is derived 
using the methods outlined above, the spontaneous breaking of chiral symmetry is accompanied 
by the appearance of a pseudoscalar massless bound state (the Goldstone boson). More details 
can be found in [75]. Consider that under a global SU(2) chiral transformation 

'f/ = e il5Tae ^ , G 1 = e G5T a e Ge n 6 r^e ? ( 31 ) 


the 2PI effective action (we ignore the dependence on the fields, as it plays no role in this 
discussion) is invariant r[G 7 ] = T[G]. One then has, for infinitesimal transformations 


W[B} = 


96T 


{ij 5 T a ,G {x',y’)} a i b i = 0 . 


SG a 'v (x',y') 

Taking one further derivative with respect to G yields 

5 T a , G (a/, y') } a /y + — ^ ~^ T a 75, cb + r a 75 , a 


(32) 


sr 


5G ab (x, y ) G a 'b> « y') 


5G ac (x, y) 


'5G cb (x,y) 


= 0 . 


(33) 


When we set G to the solution of the stationary condition, G = G^°\ the equation simplifies 


6 2 T 

$G a }) (x, 7 /) G a 'u (x 7 , y r ^j 


{ 75 r“,G (0) {x',y')} a i V = 0 . 

G=G(°) 


(34) 


If chiral symmetry is spontaneously broken then { 75 T a , G^ 0 ^ (x, y)} is non-vanishing, which means 
that there exists a solution of the BSE with pseudoscalar quantum numbers. Moreover, since 
G(°)(x,y) = G(°)(x — 7 /), the solution corresponds to one that has vanishing total momentum 
after transforming to momentum space. 

The key observation here is that, even after the effective action is truncated to some loop 
order, as long as the truncated action is invariant and both the quark self-energy and the meson 
BSE kernel are derived by taking functional derivatives of it, Goldstone’s theorem will hold 
without fine tuning. We finally note that the arguments given here translate identically to any 
global symmetry of the effective action. 


3.3. Tri-local: Three-body bound-states 

The extension of the above formulae to the three-body bound state case is rather straightforward 
[72] if one adds to the partition function a source term for the trilocal operators of interest, 
Rr S t( x ,y,z)Vr St {x,y,z). After a Legendre transformation, the effective action acquires an 
explicit dependence on the three-body vertex V. Note that depending on whether source terms 
for the proper tri-vertices of the theory are added or not, one is dealing with the 2PI or the 3PI 
effective action, supplemented with an extra vertex for the three-body bound state. 










Focusing here on the case of a baryonic bound state, we introduce sources R and R for the 
operators V — r ip r (x) r i/j s (y)'ipt(z) and V = r f> r (x) r tjjs(y)'f’t(z ), respectively. Using the stability 
arguments laid above for a would-be solution U® of the stationary condition ST/SV = 0 then 
lead to the following three-body bound state equation (see [72] for a detailed derivation) 


/ 


d A x'd A y'd A z' -^= 


<5 2 r[®,G, V] 


6V rs t{x, y, z)5V r > s ’ t ’(x', y’, z') 


G=G(°hv=v(°) 


A V r ' s i t i(x',y',z!) = 0. 


(35) 


Note that only mixed derivatives with respect to V and V do not vanish identically when one 
sets V = V(°> and V = F (0) . 

It is important to comment here that for the important cases of mixed quark-flavour states, 
and indeed mixed states in general, the procedure just outlined proceeds identically, only with 
the introduction of mixed propagators and vertices. Generally speaking, one could say that for 
each bound state of interest, one inserts the appropriate vertex in the effective action, takes 
functional derivatives with respect to it, and sets it to its vacuum value at the end. The fact 
that for flavour-diagonal mesons this vertex coincides with the quark propagator is therefore 
merely accidental. 


4. Internal structure: form factors 

The prototypical experiment for probing the internal structure of hadrons consists of (in)elastic 
scattering of a particle on a hadron. Such interactions are mediated by gauge fields, thus the 
theoretical study of hadrons must couple gauge fields so that relevant symmetries are preserved, 
as well as maintaining internal consistency of the theoretical framework. At the level of Green’s 
functions, a method coined gauging of equations exists that ensures current conservation [76-78]. 

The method of gauging can be best understood through the introduction of a new source 
term in the partition function, —J^ , {x)A^ i {x) 1 where A^ is the external, non-dynamical gauge 
field and is the current that couples to it. For example, consider the 2n-point fermionic 
Green’s function 


G (2n) (x i . ..x n ;x\ ...x' n ) = (0|T [ip 1 (x x )... ip n (xn)^ 1 (x[)... ip n (x'J] |0). (36) 

Then, its coupling to an external gauge field A M would be 

G (2n )’ M ( xi. ..x n \x’ x ...x' n ;y) = (0| T [^{xi)... ip n (xn)^ 1 (x[)... ^ n {x' n )J iX (y)\ |0), (37) 

which is called the gauged Green’s function. It can be obtained by applying a functional 
derivative with respect to A^ to the original Green’s function (we drop from now on the 
superindices indicating the number of particles) 


= 



A =0 


(38) 


and setting the external field to zero. It hence follows the ordinary rules for derivatives when 
acting, for instance, on products of functions. The proper n-body vertex, J M , for the coupling 
of the n-body system to the field A p is defined as 


= G 0 J^Go . 


(39) 


with Go again the product of full propagators. In particular, for n = 1 this defines the proper 
(fermion-photon) vertex V^ 


s» = SV^S 


(40) 







Figure 1. Gauging of the three-body kernel, as illustrated in (49). Note the minus sign 
which prevents the over counting of diagrams. 


as a result of gauging the full propagator G^ — S. A useful relation follows from gauging 
( [SS^Y and using the identity = 0 

(S-y = -S^S^S- 1 = -1^ . (41) 

For the following it is in general more convenient to work with the amputated Green’s 
function, i.e. the scattering matrix T defined in (6). Consider then a hadron described by 
the following Dyson equation 

T — K + KGqT , (42) 


with K the interaction kernels derived using the prescriptions given above. We can gauge this 
equation, following the rules of differentiation, to obtain 


T» = K tA + K^GqT + KG^T + KG 0 T» , (43) 

which can be rewritten using (7) as 

= (1 — KGq)~ 1 + K^GqT + KGqT) 

= T {K~ l K^K~ l + Gq) T . (44) 

At the bound state poles, one can introduce a bound-state electromagnetic current J ,t in a 
similar fashion as in (8) 




Vi 


rJ M - 


'h, 


P‘j +M] Pf + Mf ’ 


(45) 


where Mjj and \f \j are the initial and final bound-state masses and amplitudes, respectively. 
From (44), (9) and (10), we arrive at 


•/" = V f (Gq + GqK^Go) Vi . 


(46) 


Many of these details can be illustrated with a three-body system. Then, Go is the product 
of three full propagators S and thus its gauged analogue generates three impulse-like diagrams 


Gq = {S x S 2 S 3 y = S»S 2 S 3 + S^Ss + S&Sg = tfS 2 S 3 + S lX %S 3 + SiS 2 x 3 ■ (47) 

Here = S^Y^S 1 and the superscripts z, / denote that the propagators are evaluated before 
and after the momentum injection from the external field. The interaction kernel K must be 
decomposed into the sum of its two- and three-particle irreducible terms 


K=Y J K m) S~ l + K [3m) . 


perm. 


(48) 







Then, using (41) its gauged version is 


K “ = E tf&PliS " 1 - E K ^l)X“ + JTfspI! • («) 

perm. perm. 

It is interesting to note how the gauging procedure has automatically introduced a relative sign 
between those terms in which the external field interacts with the spectator line and those in 
which it interacts with the irreducible kernel; this ensures the absence of the over counting of 
diagrams [77,78]. The irreducible kernels themselves must be gauged once they are expressed 
in terms of the elementary degrees of freedom. A diagrammatic representation of this equation 
is shown in Fig. 1. 

The application to two-body states is entirely analogous, with the simplification that K M in 
(46) is directly the gauged two-particle irreducible kernel, K( 2 pi) m Also, the generalisation to the 
coupling of two external fields by gauging twice has been presented in [79,80]. 


5. Application to QCD 

In this section we apply the previously introduced formalism to QCD and derive the quark self¬ 
energy, quark-gluon vertex and meson Bethe-Salpeter kernel from the truncated 2PI and 3PI 
effective actions. 



Figure 2. The terms of the 3PI effective action at 3-loop relevant to the quark, quark-gluon 
vertex and meson Bethe-Salpeter kernel. 


5.1. 3PI Effective Action 

As an example, let us derive the quark self-energy, quark-gluon vertex, and meson Bethe-Salpeter 
kernel from the 3PI effective action at three-loop order. We start with the effective action, given 
in Fig. 2, wherein we have kept only those terms relevant to the discussion at hand. Note 
also that at least a three-loop expansion is required in order to obtain, upon the inclusion of a 
baryonic vertex in the effective action, a non-trivial baryon BSE kernel. 

The quark self-energy is given by (28), while solving the stationary condition 
5T[T, G, V]/5V^ b = 0 gives the quark-gluon vertex DSE 


.5T 2 [V,G,V] 

s sc* 






(50) 


0 = 


5r 2 [^,G,F] 




g=g(°) ,v=v(°) 



(51) 















In the second line of Eq.(50) we have used the vertex DSE (51) to simplify the form of the 
self-energy contribution. Note that in doing so, the self-energy is no longer a function of E, 
but rather E® which now implicitly depends upon the quark propagators G. Thus, if the BSE 
kernel is obtained by a further functional derivative of the simplified self-energy, the functional 
dependence of E® on G must be resolved. 

Thus, we use the first line in Eq. (50) to derive the Bethe-Salpeter kernel. Taking one further 
functional derivative of the quark self-energy with respect to the quark (which diagrammatically 
amounts to cutting one quark line) yields 



where we can once again use the quark-gluon vertex DSE (51) to simplify the kernel in the last 
step. 

It is interesting to note here the appearance of a ladder exchange that features two dressed 
vertices. Additionally, at this order in the truncation the doubly-dressed gluon exchange must 
necessarily be accompanied by a crossed-ladder exchange in the BSE kernel in order to preserve 
chiral symmetry and any other global symmetries of the system. 


5.2. 2PI Effective Action 

It is enlightening to compare the 3PI effective action at 3-loop to the 2PI effective action at the 
same order. We can read this off from Fig. 2 by replacing the dressed vertices with bare ones. 
Then, the quark self-energy is 


E = 


.5T 2 [tt,G] 
* SG ab 



from which the quark-gluon vertex can be inferred 



(53) 


(54) 


In a similar fashion to the above, we can write down the corresponding Bethe-Salpeter kernel 
by taking one further functional derivative of the quark self-energy with respect to the quark 




















(resolving the quark dependence of the vertex if needs be) to find 



This truncation is essentially the one employed in Refs. [81,82]; the application to baryons is 
reported in Ref. [63]. 

Notice here that (55) is structurally quite different from (52) in that the ladder exchanges 
always contain one perturbative vertex; this could be remedied by including, for example, 4-loop 
terms in the 2PI effective action. 

All of this can be compared to the rainbow-ladder truncation which follows from the 2PI 
effective action at two-loop order. Then only the two-loop terms (with vertices bare) of Fig. 2 
are required, yielding 

(56) 



Note that only this simplistic truncation lacks the flavour dependence that the two-body kernel 
necessarily features (due to implicit and explicit dependence on the quark propagator). To make 
such a truncation viable, the bare vertices are renormalisation-group improved (i.e. dressed such 
that perturbative anomalous dimensions are recovered). This accounts for the lack of interaction 
strength provided by a single gluon-exchange, and is important for both phenomenology and 
chiral dynamics. 

6. Some technical remarks 

We comment in this last section upon some technical aspects that must be taken into account 
when attempting to solve, in practice, the DSE/BSE system. 

6.1. Covariant decomposition of amplitudes 

The quantum numbers, such a spin, parity, etc. of the bound state to be studied can be 
enforced by restricting the tensor structure of the Bethe-Salpeter amplitudes to have the correct 
symmetries for these quantum numbers. 

For the description of a two-fermion bound-state, we need to provide a covariant 
decomposition ?/r^. The spinor indices of the two fermions are cq/3, while for total spin J, X is 
a product of J Lorentz indices. For spin J — 0 it is well-known that a suitable representation 
is [83,84] 

1, V*. 7 [/i 7V 1 , 7 [/ V7V 1 , (57) 

that is often simplified through the introduction of 75. Saturating these with the quark relative 
momentum k M and imposing positive parity yields — {1, ^}. Then, the general decomposition 
for a state of zero spin is 


r (k,p) = 


t 

75 


Di A± . 


(58) 

















Here, the first term selects the overall parity of the state and A± = is 

a positive/negative energy projector that introduces the (normalised) total bound-state 
momentum P. To introduce total angular momentum J, we couple the spin-zero state T to 
the two possible angular momentum tensors and 

) r ^ p ) • ( 59 ) 


These tensors are given by the traceless part of the symmetrised J-fold tensor products of a 
transversal projector transforming like a vector [85]. 

It is convenient to introduce the transverse projector Tp 1 ' — 5— P^P V j P 2 , and to denote 
its application using subscripts as follows: kp = Tp?k v \ yf = T£ OL Tp' r y 1 '. If we define the 
symmetrised J-fold tensor products 

Qui-nj = k i^ ... k ^j} ; fw-w = ... kp , ( 60 ) 


then the angular momentum tensors and T Ml are the traceless part thereof [19,24,86]. 

To describe a three-fermion bound-state we need to provide a covariant decomposition for 
^ . The spinor indices cq /?,7 correspond to the three fermionic legs, while for total spin-fc-b 1 ^, 
I is composed of one spinor and k Lorentz indices. 

Let us begin with a spin - 1 /2 baryon where 1 — 5 carries the incoming baryon spin index. 
Saturating the matrices (57) with the two relative quark momenta k M and q M and selecting 

positive parity gives Di — < 1 , ^ T , ft T (j t k Then 




/ 1 ® 1 \ 

75 <8> 75 

7T ® 7T 

7t75 ® 7t75 / 


(A (8) Dj ) ( A± 75 (7 < 8 > A + ) 


(61) 


The left tensor product denotes the outgoing quark legs with indices cq (3 and hence warrants 
the inclusion of 7 5 (7, with C — y 4 y 2 the charge conjugation matrix. The right tensor product 
describes the outgoing quark leg, index 7 , and the incoming baryon spin-index 6] the A + here 
selects the positive energy baryon. Not all elements are linearly independent; it can be checked 
that a linearly independent subspace of 64 elements can be constructed [38,87]. 

The generalisation to a state of spin -k + 1 /2 (with k integer) is obtained by extension of this 
basis 


^ a P Ul - Uk (k,q,P) = {M^-^ ® )^\k,q,P), 

( 7t"-7tA5 
njui-.-Hk — J ‘ ‘ ‘ ^75 

I 75 ’ 

with M ln " 4 ‘ k representing all combinations of k products of 7 ^% kip and qp. Here, p/u —— v k 
is the generalised Rarita-Schwinger projector. A linearly independent subspace spans 64(fc + 1 ) 
elements. As an example, we give the Rarita-Schwinger projector for spin - 3 /2 

P^i = _ pp^p , (64) 

where note that we omitted the here redundant positive energy projector A+. In this case 128 
linearly independent basis-elements can be constructed [42]. 


(62) 

(63) 





Figure 3. A sketch of the bounded parabolic region of the complex plane probed by the 
quark propagators in the Bethe-Salpeter equation. Crosses symbolise the appearance of complex 
conjugate poles in the timelike complex region. 


6.2. Euclidean spacetime and quarks in the complex plane 

Calculations using DSEs and BSEs are mainly performed in Euclidean momentum space (for 
studies using Minkowski spacetime see [88,89] and references therein). In particular, this means 
that for a bound state of mass M, the total momentum P in the rest frame is 

P= (0, 0,0, i M) . (65) 

Let us focus, for simplicity, on the case of a mesonic bound state. If the relative momentum 
between the constituent quarks is £;, the momentum of the constituents can be written as 

P± = k ± ^ • (66) 

The propagators for these constituents are functions of p±. Using (66) and assuming that the 
relative momentum k is real, one can see that 

p 2 ± — (t d= i M ) 2 (67) 

for some real and positive parameter t. This is the parametric equation for a parabola. Therefore, 
the quark dressing functions need to be known in a parabolic region of the complex plane 
(see Fig. 3). It is a general feature of the analytic structure of the quark propagators in the 
complex plane to have complex conjugate poles in the timelike complex region. These poles 
pose a limitation on the maximum mass of the bound state that can be calculated as well as 
the spacelike region for which form factors can be studied (see, e.g. [90]). The possibility of 
parametrising the quark propagators by analytic functions that allow better control over its 
singularities has been explored, e.g. in [48,91,92]. 

7. Conclusions 

In these proceedings we have collected the key aspects for a systematic study of hadronic 
properties such as masses and form-factors using the combination of Dyson-Schwinger and Bethe- 
Salpeter equations. The use of nPI effective action techniques proves to be a powerful resource 
in this respect. 

It should be clear from the presentation that the now ubiquitous rainbow-ladder truncation 
is the first term in a systematic expansion of the Bethe-Salpeter kernel using the effective action; 
in particular it appears when a renormalisation-group improved (RGI) 2PI effective action at 
two-loops is used. Present investigations [63] are aimed at extending the calculation of meson 
and baryon spectra using kernels derived from a three-loop truncation of the effective action. 
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